Quantum Phase Transitions in the Ising model in spatially modulated field 
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The phase transitions in the transverse field Ising model in a competing spatially modulated 
(periodic and oscillatory) longitudinal field are studied numerically. There is a multiphase point 
in absence of the transverse field where the degeneracy for a longitudinal field of wavelength A is 

(1W5)2JV/A for a system with N 

spins, an exact result obtained from the known result for A = 2. 
The phase transitions in the V (transverse field) versus ho (amplitude of the longitudinal field) 
phase diagram are obtained from the vanishing of the mass gap A. We find that for all the phase 
transition points obtained in this way, A shows finite size scaling behaviour signifying a continuous 
phase transition everywhere. The values of the critical exponents show that the model belongs to 
the universality class of the two dimensional Ising model. The longitudinal field is found to have 
the same scaling behaviour as that of the transverse field, which seems to be a unique feature for 
the competing field. The phase boundaries for two different wavelengths of the modulated field are 
obtained. Close to the multiphase point at h c , the phase boundary behaves as (h c — ha) h ', where b 
is also A dependent. 
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PACS numbers: 64.60.Cn, 64.60.Fr, 42.50.Lc 



I. INTRODUCTION 



Competitive interactions in magnetic models can gen- 
erate frustration leading to the existence of a large num- 
ber of degenerate states Q. For example, in the ANNNI 
(axial next nearest neigbour Ising) model [jl|-§| , the pres- 
ence of a second neighbour antiferromagnetic interaction 
induces a competition. In absence of any fluctuation, 
a large number of degenerate states exist at a so called 
multiphase point. In spin glasses interactions are 
random, and the mixture of ferromagnetic and antiferro- 
magnetic interactions result in frustration and huge de- 
generacy. 

Phase transitions, both classical and quantum, in 
many of these systems are well studied jL)-f§. While the 
short range Ising like models do not have any finite tem- 
perature phase transition in one dimension, the phase 
diagram in higher dimensions often shows very rich struc- 
tures. The quantum phase transitions are present even 
in one dimension. It is interesting to observe that the 
zeroth order thermal fluctuation in two dimensions and 
correspondingly zeroth order quantum fluctuations in one 
dimension can destroy all order at the fully frustrated 
multiphase point (|J|] in ANNNI-like models. 

The Ising models in transverse field form a special class 
of quantum systems and to the best of our knowledge, 
have different features in comparison to quantum spin 
models in which the quantum effect comes through the 
cooperative interaction. 

In this paper, we will discuss a quantum Ising model 
with competition, where the competition is generated by 
the presence of an external longitudinal field which is spa- 
tially periodic. We choose the periodicity of the field to 
be commensurate with the lattice periodicity. The field is 
also oscillatory so that the net field on the entire system 



sums out to zero. Classically, the ferromagnetic Ising sys- 
tem will not undergo any phase transition in a uniform 
field. On the other hand, spatially modulated fields can 
play important roles in critical phenomena. Classical sys- 
tems in periodic potential (Frenkel-Kontorova model), for 
example, show very interesting features leading to novel 
concepts like the devil's staircase |fj . 

Some interesting and relevant features of the model in 
the classical limit are discussed in section II. The phase 
transitions in the ferromagnetic model in the modulated 
field in one dimension at zero temperature driven by 
quantum fluctuations generated by the transverse field 
are described in section III. The results are discussed and 
concluding remarks are made in section IV. 



II. THE MODEL: CLASSICAL GROUND STATE 
AND DEGENERACY 
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The Hamiltonian for the system with N spins is given 
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(i) 



The form of the longitudinal periodic field hi is like 
this: hi = ho at i = (nX + 1) to (2n + l)A/2 (n — 
0, 1, N/X — 1) and — ho elsewhere. The wavelength 
of the field is denoted by A. T is the strength of the 
transverse field. 

Note that the modulated field can be chosen to be of 
various forms, for simplicity we will choose a square wave 
form. While the ferromagnetic interaction tries to align 
the spins parallely, the periodic field will try to modulate 
the system spatially. We will consider periodic boundary 
condition. 
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We will first discuss some known results in the con- 
text of our model in the classical limit (i.e., T = 0). The 
classical Hamiltonian is 

N N 

H d = -J*^2 - X! h^i- ( 2 ) 

i=l i=l 

Let us first take the case A = 2: the field is positive on 
odd sites and negative on even. Now one can simply use 
a transformation — * —£>2i- This effectively makes the 
Hamiltonian 

N N 

H c i = SiSf+i ~ Sf, (3) 

i=l i=l 

i.e, an antiferromagnetic system in a uniform field |||]] 
which is solved exactly by transfer matrix method. 

Again, if we consider the ANNNI chain, the Hamilto- 
nian is: 

JV N 

H ANNNI = ~Jl SfSf +1 + S?Si+2- ( 4 ) 

i=l i=l 

The above can be cast to the form of (3) by the transfor- 
mation SfS* +l — > S* [Q with Ji playing the role of ho 
and J2 the role of J. Note that J2 is positive, describing 
antiferromagnetic interaction by definition. The sign of 
Ji is unimportant. 

The two models described by equations (3) and (4) 
are well studied |^|,0J3] and they have similar classical 
ground states: upto J /ho = 0.5 (or J2/J1 = 0.5), the sys- 
tem is antiferromagnetic (ferromagnetic for ANNNI when 
J\ > 0), and paramagnetic (modulated for ANNNI) be- 
yond this value. An infinite number of degenerate states 
exist at this point. The degeneracy can be calculated ex- 
actly [HJ|] and is equal to t n where r = (1 + v / 5)/2. In 
the language of the original classical model defined in (2) 
with A = 2, one has ferromagnetic phase upto ho/ J = 2 
and paramagnetic phase for ho/ J > 2. 

Even when A 7^ 2, one can rescale the system by a fac- 
tor of A/2 so that the Hamiltonian is effectively written 
as in (3). It can easily be checked that the multiphase 
point is now at h c / J — 4/A and the degeneracy given by 

T 2N ' X . 

It is extremely important to note here that although 
the antiferromagnetic and ferromagnetic systems have 
identical critical behaviour, the presence of a field makes 
things drastically different. Therefore the modulated 
field in the ferromagnetic system is effectively a uniform 
field in the mapped antiferromagnetic system. We will 
specifically take the interaction to be ferromagnetic in 
order to avoid any ambiguity or confusion. 

III. QUANTUM PHASE TRANSITION AND 
PHASE DIAGRAM AT T=0 

There are two limits in which the results are exactly 
known for the quantum Hamiltonian (1) at zero temper- 



ature. The classical limit r = is already discussed. 
In the limit ho = 0, the problem is also exactly solved 
JlOj l : there is a phase transition at a finite value of the 
transverse field: T c (0)/J = 1. This is a quantum critical 
point and the phase transition is continuous in nature. 
In a spatially periodic and oscillatory longitudinal field 
which is competing in nature, the phase transition should 
occur even at lower values of the transverse field, T c (ho). 

Vanishing of the mass gap A(= E\ — Eq where E\ and 
Eq are the energies of the first excited and the ground 
state respectively) in the thermodynamic limit is a sig- 
nature of a continuous phase transition. We obtain the 
mass gaps for finite chains of length L (N = L d = L, 
where d the dimensionality is one here). 

We employ Lanczos' method to diagonalise the Hamil- 
tonian matrix which is calculated in the basis in which 
S z is diagonal. The field being spatially modulated and 
also commensurate with the lattice periodicity, the sys- 
tem sizes are restricted to specific values for a given wave- 
length. On the other hand, one requires data for different 
system sizes for a finite size scaling procedure. Since it 
is difficult to work with large sizes in a diagonalisation 
scheme anyway, we have obtained the phase diagrams for 
A = 2 and 4 only. The maximum size taken for both cases 
is L = 16. The result for a very large value of A can be 
easily guessed. 

We indeed find that A vanishes with 1/L not only 
for small values of ho but right upto the critical field 
h c = 4J/A. The critical points r c (h ) can be obtained 
by extrapolation, and a more accurate estimation is done 
using finite size scaling of the mass gap, which shows 
scaling behaviour. We defer a detailed discussion on the 
scaling of mass gaps and nature of transitions and first 
summarise the results for the phase diagram. 

In Fig. 1, the phase diagrams for A = 2 and 4 are 
shown. We have calculated the ferromagnetic order: 
(to) = j i '}2 i {S z ) as well as the "staggered magnetisa- 
tion" defined as 4 XX^? hi) /ho- The phase boundary 
we obtain separates a ferromagnetic phase with nonzero 
(to) and a paramagnetic phase with (m) = 0. However, 
since the longitudinal field is always present, the stag- 
gered magnetisation is nonzero all over the paramagnetic 
phase. 

It is also observed that as soon as h c /J = 4/A is 
reached, the system behaves as if only a field is present 
and therefore no phase transition is observed here for any 
non-zero T: the mass gap remains finite in the thermody- 
namic limit almost without any system size dependence. 
As in some other models with competing interactions, 
in this model also, the zeroth order quantum fluctuation 
destroys the order at h c and lifts the degeneracy 

The entire phase boundary could not be fit to any sim- 
ple form; we could obtain the behaviour close to the 
multiphase point quite accurately for both A = 2 and 
A = 4. Here the phase boundary appears to be of the 
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form (4/A - h /J) b where b = 0.77 ± 0.02 for A = 2 and 
b = 0.49 ± 0.01 for A = 4. 

For A — > oo, the critical field h c — > 0, the ferromagnetic 
phase exists only for ho = 0, and the phase boundary will 
coincide with the ho = axis. The exponent b decreases 
with A, ultimately vanishing for A — > oo. 

We end this section with a discussion on the scaling 
behaviour of mass gaps and the nature of phase transi- 
tions. In zero longitudinal field, the finite size scaling for 
the mass gap gives: 

A ~ L-*/((T - T^L 1 '") (5) 

where z is the dynamic exponent and v is the correlation 
length exponent. 2 = 1 and v = 1 for the one dimen- 
sional transverse field Ising model. The latter can be 
mapped to a two dimensional classical Ising model. / is 
a universal scaling function. 

We find the above finite size scaling form to be valid for 
non-zero values of ho in the sense that when LA is plot- 
ted against (T — T c (ho))L l / v the data for different system 
sizes again collapse for any value of ho (see Fig 2). The 
data collapses are obtained with z = 1 and v = 1 for 
all ho, indicating that the model belongs to the univer- 
sality class of the two dimensional classical Ising model 
even for non-zero longitudinal modulated fields. As ho is 
increased, the collapse gets limited to a smaller region. 
The scaling function depends on the value of ho- 

The phase transition points depend on the values of 
both r and ho- The plots in Fig. 2 are done by keeping 
ho fixed and varying T around the critical value. One 
can study the scaling the other way also: keep T fixed 
and vary ho- We find that there is a collapse of data for 
different sizes again, as shown in Fig. 3. The striking 
feature is that the scaling argument is (ho — h c {T))L 1 / v 
with v = 1. This shows that the competing longitudinal 
field scales in the same way as does the transverse field. 
The results for A = 4 show the same scaling behaviour as 
in Figs 2 and 3. We discuss this feature again in section 
IV. 

One may conclude from the above observations that 
the phase transitions are continuous everywhere (except 
at r = 0) as the mass gaps vanish continuously and show 
scaling behaviour as in conventional critical phenomenon. 



IV. DISCUSSIONS AND CONCLUSIONS 

The Ising model in the presence of a transverse field 
and a longitudinal field cannot be solved exactly. Re- 
cently, most of the studies on transverse Ising models 
have involved randomness in the transverse field, inter- 
actions or longitudinal field. We have considered a very 
simple model in which there are both quantum effects 
and frustration, where the frustration comes as a result 
of the modulated nature of the external field and no ran- 
domness is involved. 



In absence of the transverse field, the system has a mul- 
tiphase point. The degeneracy at the multiphase point 
in the absence of fluctuations is easily calculated for any 
A from the known result corresponding to A = 2. 

We have obtained the phase diagram of the model in 
the r — ho plane at zero temperature using numerical 
methods. Our results show that a continuous phase tran- 
sition occurs everywhere except at the multiphase point 
where r c = and a first order phase transition is known 
to exist. The transition at this point is of first order in the 
sense that the order parameter (magnetisation) discon- 
tinuously vanishes at this point. The values of the critical 
exponent z and v obtained in this model are identical to 
those of the transverse Ising model, showing that the pe- 
riodic longitudinal field is irrelevant. Thus it belongs to 
the classical two-dimensional Ising universality class. 

Our conclusions are based on the scaling behaviour of 
the mass gap. As in conventional critical phenomena, it 
shows finite size scaling behaviour close to the critical 
point. It is to be noted that according to scaling argu- 
ments in classical critical phenomena, there are only two 
relevant fields for magnetic systems, temperature (T) and 
magnetic field (h), and the correlation length scales as: 

£ = L g{(L l ' v t > hL^ + ^' v ). (6) 

The deviation from the critical temperature is denoted 
by t, [3 and 7 are the critical exponents associated with 
the order parameter and susceptibility respectively and 
g is a universal scaling function. Here the field h and the 
temperature have different scaling dimensions. There is, 
however, no criticality associated with the uniform or- 
dering field h. In the present model, the two important 
quantities are the transverse field and the competing lon- 
gitudinal field. The quantum fluctuations through T may 
be compared to the thermal fluctuations. The roles of T 
and T may be considered to be equivalent; however, the 
roles of the competing field and ordering field are com- 
pletely different. That is why a scaling form as in (^) will 
not be valid here. On the other hand, from the study of 
the finite size scaling of the mass gaps (section III), it 
appears that the competing field and the transverse field 
have the same scaling behaviour. We claim this to be a 
unique feature of the non-ordering competing field. 

Beyond the multiphase point, the field is dominat- 
ing, and no continuous phase transition can exist here; 
the mass gap remains finite in the thermodynamic limit. 
This is in contrast to the competitive models like ANNNI 
model, where one obtains phase transitions beyond the 
fully frustrated point as well. This difference is of course 
due to the different roles played by interaction and field. 
The correspondence between the Ising model in compet- 
ing field and the ANNNI model mentioned in section II 
is valid only in one dimension and in the classical limit. 

From the present study, however, existence of other 
modulated phases, separated by first order lines, cannot 
be ruled out. Such a possibility is strongest close to the 
multiphase point. A closer inspection, eg., using pertur- 
bation methods could be done from which it may also be 
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possible to verify the behaviour of the phase boundary 
and get an estimate of b defined in section III. 

A few other points need to be mentioned. In the ran- 
dom field transverse model, a tricritical point was ob- 
tained |ll]] . The phase diagram of certain antiferromag- 
nets like FeCl2 in a uniform field also show a tricritical 
point. In comparison, our model does not have any such 
point. Rather, the phase boundary to the left of the mul- 
tiphase point is very similar to that of the ANNNI model 
in higher dimensions. In the latter also, the next neigh- 
bour interaction cannot change the universality class of 
the Ising model. The mapping to the antiferromagnetic 
model in the classical limit and subsequent description 
of the ground state assumes that we are considering hy- 
percubic lattices only, where, in absence of any field at 
zero temperature, the antiferromagnet will order. This is 
not true in lattices in which the antiferromagnetic system 
is geometrically frustrated (e.g., in a triangular lattice); 
such systems in transverse (and longitudinal) field in two 
dimensions have been recently studied [Q . 

We have considered only a square wave form for the 
modulated field in this discrete model. The features 
should not be different for a different form, however, 
there will be quantitative changes. In a continuum 
model, it may be more convenient to choose a sinu- 
soidal form. More complex spin patterns are obtained in 
(antiferromagnetic) systems with long range interactions 
in a field, e.g ., the long range antiferromagnet in a uni- 
form field |L3|] exhibits a complete devil's staircase. It is 
expected that there will also be very interesting features 
if the modulation in the field is incommensurate with the 
lattice periodicity. 
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FIG. 1. The phase boundaries between ferromagnetic and 
paramegnetic regions are shown for the transverse Ising model 
in a spatially modulated field for two wavelengths A = 2 
and 4. Near h c (X)/J = 4/(A), the boundaries fit to the 
form (h c (\) — ho) b , shown by the continuous curves, where 
b = 0.77 ± 0.02 and 0.49 ± 0.01 for A = 2 and A = 4 respec- 
tively. 

FIG. 2. The scaling plot for the scaled mass gaps L z A with 
(r — Tciho))!}/" for A = 2. The four sets correspond to the 
values of ho = 1.9, 1.8, 1.6 and 1.0 from top to bottom, each 
showing the collapse of data points for L = 8, 10, 12 and 16 
with z = 1 and v = 1. 

FIG. 3. The scaling plot for the scaled mass gaps L z A with 
(ho — /i c (F))L 1 ^ for A = 2. The four sets correspond to the 
values of V — 0.25, 0.45, 0.65 and 0.95 from top to bottom, 
each showing the collapse of data points for L — 8, 10, 12 and 
16 with z — 1 and v = 1. 
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Fig 1: Phase diagram 




0.5 1 1.5 2 2.5 

ho/ J 



Fig2: Scaling with T — r c 
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Fig3: scaling with (h — h c ) 
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